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The influence of a local anisotropy of random orientation on a ferromagnetic phase 
transition is studied for two cases of anisotropy axis distribution. To this end a 
model of a random anisotropy magnet is analyzed by means of the field theoretical 
renormalization group approach in two loop approximation refined by a resum- 
mation of the asymptotic series. The one-loop result of Aharony indicating the 
absence of a second-order phase transition for an isotropic distribution of random 
anisotropy axis at space dimension d < 4 is corroborated. For a cubic distribution 
the accessible stable fixed point leads to disordered Ising-like critical exponents. 



1 Introduction 

Modern understanding of universal properties of matter in the vicinity of crit- 
ical paints is mainly due to the application of the renormalization group (RG) 
ideasJd Applied to the problems of condensed matter physics in the early 1970s, 
the RG technique proved to be a powerful tool to study critical phenomena. 
For example, expressions for critical exponents governing the magnetic phase 
transition in regular systems are known bu_now with record accuracy both for 
isotropico (O(m) symmetrical) and cubida magnets. The RG approach also 
sheds light on the influence of structural disorder on ferromagnetism. In the 
present article we will apply the field theoretical RG approach to study pe- 
culiarities of magnetic behaviour influenced by disorder in a form of random 
anisotropy axis.u It is special pleasure for us to dedicate this paper to Prof. 
Hagen Kleinert on the occasion of his 60th anniversary. His contribution to 
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the field is hard to be overestimated. 

Although an influence of a weak quenched structural disorder on universal 
properties of a ferromagnetic phase transition has already been a problem 
of intensive study for several decades, there remains a number of unsettled 
questions. Here, one should distinguish between random site, random-field 
and random anisotropy magnets. A weak quenched disorder preserves second- 
order phase transition in three-dimensional (d = 3) xandom site magnetsQ but 
can destroy this transition in random field systemsO for d < 4. The situation 
for the random-anisotropy magnets is not so clear. 

Typical examples of random-anisotropy magnets are amorphous rare- 
earth - transition metal alloys. Some of these systems order magnetically 
and for the description of the ordered structure it has been proposeda to 
consider a regular lattice of magnetic ions, each of them being subjected to 
a local anisotropy of random orientation. The Hamiltonian of this random 
anisotropy model (RAM) readsH 

H = - Jr,r'SrSW - A> J](£rSr) 2 , (1) 

where Sn is an m-component vector on a lattice site R, Jr.r' is an exchange 
interaction, D is an anisotropy strength, and % is a unit vector pointing in 
the local (quenched) random direction of an uniaxial anisotropy. 

The model has been investigated by a variety of techniques including 
mean-field theoryjj comput er simulations,!] l/m-expansion,Q renormalization 
group e-expansion. The talit. case of an infinite anisotropy has been 

subject of a detailed study as well.lialiJ However the nature of low-temperature 
phase in RAM is not completely clear up to now, althoughj-scvcral low- 
tempeiature phases were discussed like ferromagnetic orderingjZHa spin-glass 
phasep'Q and quasi long-range ordering.113 

The nature of ordering is connected with the distribution of the ran- 
dom variables ir in (Q). For an isotropic distribution arguments similar to 
those applied by Imry and Mali3 for a random-field Ising modeLJjH±ag about 
the absence of ferromagnetic order for space dimensions d ^4,E3li3 whereas 
anisotropic distributions may lead to a ferromagnetic order x3 

Application of Wilson RG technique to RAM with the isotropic distribu- 
tion of a local anisotropy axis suggestalJ the possibility of "runaway" solutions 
of the recursion equations. Such a behaviour was interpreted as a smeared 
transition. However this result was obtained in first order of the e-expansion 
and remains to be confirmed also in higher orders. 

Here, we will report results obtained by means of the field theoretical RG 
technique in two loop approximation refined by a resummation of the resulting 
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asymptotic series. We will consider two cases of distribution of the random 
anisotropy axis and show that a ferromagnetic second-order phase transition 
takes place only when the distribution is non-isotropic. Moreover we will 
show that the RAM provides another example of a disordered model, where 
the only possible new critical behaviour is of "random Ising" type, similar. to 
the site diluted magnets. □ More detailed results can be found in Refs£3'E!l 



2 Isotropic case 

In order to deal with quenched disorder, one way to obtain the effective Hamil- 
tonian of a RAM is to make use of the replica trick. For a given configuration 
of quenched random variables ccr in (Q) the partition function may then be 
written in a form of functional integral of a Gibbs distribution depending on 
xr. To average over configurations one should complete the model by choos- 
ing certain distribution of £r. We will analyze two cases: first the isotropic 
case, where the random vector x points with equal probability in any direc- 
tion of the m-dimensional hyperspace, and second the cubic case, where x lies 
along the edges of m-dimensional hypercube (cubic case) . Other distributions 
may be considered as well. In the first case the distribution function reads: 

*>-(/«)"'- w 

Following the above described program, one ends up with the replica n — > 
limit of the effective Hamiltonianty: 

n eff = -J d d R\\[^\ 2 +\v0\ 2 ] + U0 M 4 +^|<n 4 + 



n 



<<•(> Y. »>;o/o/>. en 

where fig 2 and uq, vq, u)q are defined by Dq and familiar bare couplings of an 
m-vector model, and 4> a = </> R is a m-dimensional vector, \<p\ 2 — y~) ^ |^ Q | 2 . 
The bare couplings are restricted by uq > 0, vq > 0, wq < 0. Furthermore, 
values of uq and wq are related to appropriate cumulants of the distribution 
function (^) and their ratio equals wq/uq = — m. Note that the symmetry of 
wo and Vo terms corresponds to the random site m-vector model.Ell However 
the uo-term has an opposite sign. 

In order to study long-distanca properties of the Hamiltonian (|3[), we use 
the field-theoretical RG approach.B Here the critical point of a system corre- 
sponds to a stable fixed point (FP) of the RG transformation. We apply the 
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massive field theory renormalization schemeE2l performing renormalization at 
fixed space .dimension d and zero external momenta. In two-loop approxima- 
tion we getlij expressions for the RG functions in form of asymtotic series in 
rcnormalizcd couplings u,v,w. 
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Figure 1. Fixed points of the RAM with isotropic distribution of a local anisotropy axis. 
The fixed points located in the octant « > 0, v > 0, w < are shown. The filled box shows 
the stable fixed point, the cross denotes typical initial values of couplings. 



As it was mentioned in the introduction, the only known RG results for 
RAM with isotropic distribution of the local anisotropy axis so far are those 
obtained in first order in e.l!3 In total one gets eight fixed points. All FPs with 
u> 0, v > 0, w <0 appear to be unstable for e > except of the "polymer" 
0(n — 0) FP III which is stable for all m (see Fig. Q). However the presence 
of a stable FP is not a sufficient condition for a second-order phase transition. 
In order to be physically relevant, the FP should be accessible from the initial 
values of couplings. This is not the case for the location of FPs shown in Fig. 
[l]. Indeed starting from the region of physical initial conditions (denoted by 
the cross in the Fig. [I]) in the plane of v = one would have to cross the 
separatrix joining the unstable FPs I and VI. This is not possible and so one 
never reaches the stable FP III. As far as both FPs I and VI are strongly un- 
stable with respect to v, FP III is not accessible for arbitrary positive v either. 
Finally, runaway solutions of the RG equations show that the second-order 
phase transition is absent in the model. The main question of interest here 
is whether the above described picture of runaway solutions is not an artifact 
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Table 1. Resummed values of the fixed points and critical exponents for isotropic case 
in two-loop approximation for d = 3. We absorb the value of a one-loop integral into the 
normalization of couplings. 
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of the e-expansion. To check this, we used a more refined analysis of the FEs 
and their stability, considering the series for RG functions directly at d — 3E3 
It is known that series of this type are at best asymptotic and a resummation 
procedure has to be applied to obtain reliable data on their basis. We make 
use of Pade-BoraLresummation techniques ,E3 first writing the RG functions 
as resovent seriesnj in one auxiliary variable and then performing resumma- 
tion. Numerical values of the FPs are given in Table [j]. Resummed two-loop 
results qualitatively confirm the picture obtained in first-order e-expansion: 
stability of the FPs does not change after the resummation. This supports the 
conjecture of AharonyEj that accessible stable FP for the RAM with isotropic 
distribution of the local anisotropy axis is absent. In the table, we list values 
of the correlation length and pair correlation function critical exponents v and 
r\ which are resummed in a similar way. As they are calculated in unstable 
FPs they have rather to be considered as effective ones. 



3 Cubic case 



Let us now consider the second example of anisotropy axis distribution, when 
the vector £r (Q) points only along one of the 2m directions of axes ki of a 



ram: submitted to World Scientific on February 1, 2008 5 



cubic lattice : 

1 m 

P(x) = —J2[ s{m) (x-h) + SW(x + k)}. (4) 



2m 

i=l 



The rationale for such a choice is to mimic the situation when an amorphous 
magnet still "remembers" the initial (cubic) lattice structure. Repeating the 
procedure described in the previous section, one ends up with the following 
effective Hamiltonian which is of interest in the limit n — > 0E9: 



H 



eff 




i=l a,/3=l i=l a=l J 

Here, the bare couplings are ito > 0, «o > 0, wq < 0. The yo term is gener- 
ated when the RG transformation is applied and may be of either sign. The 
symmetry of wq terms differs in (||) and (||). Furthermore, values of wq and 
i*o differ for Hamiltonians (||) and (^) but their ratio equals — m. again. 

We apply the massive field theory renormalization schemec 2 ] and get the 
RG functions in two-loop apprpximation.tj As in the previous case we repro- 
duce the first-order e-results.E3 Now one gets 14 FPs. However, in first order 
of the e-expansion all FPs with u > 0,v > 0,w < appear to be unstable 
for e > except of the "polymer" 0(n = 0) FP HI which is stable for all m 
but not accessible (see Fig. ||) . Now the account of the £ 2 -terms qualitatively 
changes the picture. Indeed, the system of equations for the FPs appears 
to be degenerated at the one loop level. As known from other cases in two 
loop order this leads to the appearance of a new FP which is stable and is 
expressed by a-^/e series£d The possibility of such scenario was predicted al- 
ready in Refs.Ea However it remained unclear whether there exist any other 
accessible stable FPs. 

Applying Pade-Borel resummation we get 16 FPs. Values of the FPs 
with u* > 0, v* > 0, w* < are listed in the Table |[ The last FP XV in 
Table |^ corresponds to the stable FP of -^/e-expansion. It has coordinates 
with u* — v* — 0, w* < and y* > and is accessible from the typical 
initial values of couplings (shown by a cross in the Fig. ||). Applying the 
resummation procedure, we have not found any other stable FPs in the region 
of interest. The effective Hamiltonian (]^) at u = v = in the replica limit 
n — > reduces to a product of m effective Hamiltonians of a weakly diluted 
quenched random site Ising model. This means that for any value of m > 1 
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Figure 2. Fixed points of the RAM with distribution of a local anisotropy axis along 
hypercube axes for v = 0. The only fixed points located in the region u > 0, w < are 
shown. Filled boxes show the stable fixed points, the cross denotes typical initial values of 
couplings. 



the system is characterized by the same set of critical exponents as those of 
a weakly diluted quenched random site Ising model. In the Table, we give ^| 
values of critical exponents in the other FPs as well: if the flows from initial 
values of couplings pass near the these FPs, one may observe an effective 
critical behaviour governed by these critical exponents. 

4 Conclusions 

We applied the field theoretical RG approach to analyze critical behaviour of 
a model of random anisotropy magnets with isotropic and cubic distributions 
of a local anisotropy axis. The origin of a low-temperature phase in this 
model is not completely clear. General arguments based on an estimate of 
the energy for formation of magnetic domainal^ lead to the conclusion that for 
d < 4 a ferromagnetic order is absent. E3'EZl However, these arguments do npt 
take into account entropy which may be important for disordered systemsO 
Furthermore, they do not apply for anisotropic distributions of the random 
axis £3 

In the RG analysis the absence of a ferromagnetic second-order phase 
transition corresponds to the lack of a stable FP of the RG transformation. 
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Table 2. Resummed values of the FPs and critical exponents for cubic distribution in 
two-loop approximation for d = 3. We absorb the value of a one-loop integral into the 
normalization of couplings. 
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However in the case of RAM with isotropic distribution of a local anisotropy 
axis the scenario differs. Our two-loop calculation leads to a 0(n = 0) sym- 
metric FP which is stable for any value of m for both isotropic and cubic 
distributions of a random anisotropy axis. Note that this FP is not accessible 
from the initial values of the couplings. We checked the location of the FPs up 
to second order in the e-expansion and by means of a fixed d = 3 technique 
refined by Pade-Borel resummation. 

In the case of isotropic distribution aLa random anisotropy axis our anal- 
ysis supports the conjecture of Aharonyld based on results linear in e about 
runaway solutions of the RG equations. For the cubic distribution we get two 
stable FPs. One of them (FP III in Fig. |^) is not accessible as in the isotropic 
case. But the disordered Ising-like FP (FP XV in Fig. |^) may be reached 
from the initial values of couplings. Applying the resummation procedure we 
have not found any other stable FPs in the region of interest. This means that 
RAM with cubic distributions of a random anisotropy axis is governecLhy. a 
set of critical exponents of a weakly site diluted quenched Ising model.ErEj 

To conclude we want to attract attention to a certain similarity in the 
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critical behaviour of both random-siteEl and random-anisotropyu quenched 
magnets: if there appears a new critical behaviour it always is governed by 
critical exponents of site- diluted Ising type. The above calculations of a "phase 
diagram" of RAM are based on two-loop expansions improved by a resumma- 
tion technique. Once the qualitative picture becomes clear there is no need to 
go into higher orders of a perturbation theory as far as the critical-exponents 
of site-diluted Ising model are known by now with high accuracy.E3 
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